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"^Abstract 
C . 

Oh' Using an extensive compilation of literature supernova rate data we study to which extent its evolution constrains the star 
^ formation history, the distribution of the type la supernova (SNIa) progenitor's lifetime, the mass range of core-collapse supernova 
5^ (CCSN) progenitors, and the evolution of the iron density in the field. 

^ ' We find that the diagnostic power of the cosmic SNIa rate on their progenitor model is relatively weak. More promising is the 
Vise of the evolution of the SNIa rate in galaxy clusters. We find that the CCSN rate is compatible with a Salpeter IMF, with a 
""^minimum mass for their progenitors > 10 Mq. We estimate the evolution in the field of the iron density released by SNe and find 
that in the local universe the iron abundance should be ~ 0.1 solar. We discuss the difference between this value and the iron 
abundance in clusters. 
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1. Introduction 

, As the last stage of the evolution of stars in a wide range 
of masses, supernovae play a unique role in the cycling of 
piatter in galaxy evolution. They release in the interstellar 
medium all of the iron available in the Universe, which is 
produ ced (mainly as radioa ctive Nickel 56) in the explosion 



itself (jWooslev et al.l 120021 ). The supernova rate (SNR) is 



then closely related to the rate of Fc enrichment of galaxies 
and of the universe in general. 

There are two main families of supernovae, the "Core- 
Collapse Supernovae" (CCSN) which regroup into the ob- 
servational types II, lb and Ic, and the type la SNe (SNIa), 
which are the product of the thermonuclear explosion of a 
white dwarf (WD), accreting matter from a companion. 

CCSN are the result of the death of massive stars. Mod- 
eling the explosion is rather difficult due to the complicated 
physical processes taking place i n the extreme co nditions 
of the collapsing F e core (see e.g. lWooslev fc Weavcr„1986c 
Ijanka et al. 2007 ). This reflects into an intrinsic difficulty 
to predict the exact mass range of the progenitors of the dif- 



ferent CCSN kinds, as well as of the amount of radioactive 
nickel released to the interstellar medium by each event. 
However, there is no doubt that the lifetimes of CCSN pro- 
genitors are short, less then a few 10 Myr, which implies 
that the CCSN rate closely follows the star formation rate. 

The main frame of the SNIa theoretical paradigm needs 
a CO white dwarf in a close binary system. The primary 
star results from the evolution of an intermediate mass star 
(M < 8 Mq) which develops a degenerate CO core. When, 
due to evolution, the secondary expands and fills its Roche 
Lobe, the CO WD may accrete and grow in mass, and, if 
the Chandrasekhar limit is reached, C is ignited under de- 
generate conditions, initiating a thermonuclear explosion. 
This is the Single Degenerate (SD) ch annel for SNIa's pro- 



duction (see e.g. Hachisu et al. 19991 ). Conversely, if effi- 
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cient growth of the CO WD does not occur, the envelope 
of the secondary is dispersed in the interstellar medium, 
leaving behind two close WDs. The binary system has a 
second chance to give rise to a SNIa event, since, due to the 
emission of gravitational waves radiation, orbital energy is 
lost and the two WDs eventually merge. If the total mass 
of this Double Degenerate (DD) system e xceeds the Chan- 
drase khar limit a SNIa may be produced (|lben fc Tutukov 
19841) . While the modeling of the binary system evolution 



from its birth to the final explosion is subject to many rel- 



ativcly uncertain prescriptions, the nucleosynthetic prod- 
uct of the explosion seems more robustly assessed, and it 
consists of mainly Fe peak elements, and in termediate (Si 
group) mass elements ( Iwamoto et al.l[l999f ). The time de- 
lay between the birth of the binary system and the final 
explosion spans a wide range, from a few 10^ yrs, corre- 
sponding to the nuclear evolutionary lifetime of the more 
massive systems, up to a Hubble time or more, especially 
for the DD systems born with the larger separations. As 
a consequence, the scaling of the SNIa rate with time will 
reflect not only the star formation history of the parent 
stellar system, but also the distribution of the delay times. 

In this paper we analyze the constraints on the pro- 
genitors of CCSN, on the distribution of the delay times 
of SNIa. and on the cosmic history of star formation 
which can be derived from the measurements of the 
cosmic SN rates of the two kinds. Previo u s studies on 
this s u bject include Jorgensen et ahl (1997 ): Sadat et al 



(Il998l); iDahlen fc FranssonI (l999l) : iDahlen et all (|2004 ): 



Forster et al.l (|2006l ). With respect to these, we try to test 



the galaxy population in cluster seems rather homogeneous 
and old. First measurement in high redshift clusters has 
been made b y lCal-Yam et al. ( 2002 ). with little statistics. 
ReissI (|2000h used the Mount Stromlo Ahell Cluster SN 



Search to derive a rate within clusters. Th e latest mea- 
surem ent of this kind has been published bv ISharon et al 



(|2007l ). where a summary of all current SNIa rate measure- 
ments in clusters is presented. 

First high-redshift measurements of the CCSN rate have 
been published by two teams almost si multaneously, the 
STRES S survev (ICaDDellaro et al.l l200,ql) and the GOODS 



survey (IDahlen et al 



. 20041) The result from STRESS 

survev bv iBotticella et all (|2008f) supersedes the one of 
Cappellaro et al. ( 2005f ). Tab.[2lsummarizes these measurc- 



the conclusions versus the choice of relevant parameters, 
by considering a variety of options for the distribution of 
the delay times and for the SFH. 

The paper is organized as follows: in Sec. [2l SNR mea- 
surements are reviewed; in Sec. |3l the various theoretical 
ingredients hidden in the SNR evolution are dissected; in 
Sec. S] we discuss the fit of the models to the observed cos- 
mic SNIa rate in the field and in clusters, and to the ob- 
served cosmic CCSN rate. In Sec. [5] we discuss the impli- 
cations of our results on the evolution of the Fe density in 
the universe. Conclusions are drawn in Sec. [H 

We use the concordance cosmological model (rij\/„ , I^Ao ) = 
(0.3, 0.7); with a Hubble constant set to Ho ~ 7O/170 kms~-'^Mpc 



ments. The high-redshift CCSN sample used to derive the 
rate is no t completely based on sp e ctroscopic confirmation , 
bot h forlCappellaroet"an (|2005[ ): IBotticella et all (|2008l ) 
and Dahlen et al. ( 2004 ) measurements. Another source of 
uncertainty may come from the extinction within the host, 
since CCSN are known to occur shortly after their birth in 
highly du sty star formation region s of galaxies. Measure- 
ments by ICappellaro et all (Il999h are not corrected for 



Dahlen et al. I (I2OO4I)" attempted to make such a 



extinction. 

correction by simulating the host dust distribution accord- 
ing to the di s k incl ination, following the prescription by 



Hatano et al.l (|l998l ) . Their corrected rates (as used in this 
paper) are a factor of ~ 2 higher than the uncorrected. We 
notice, however, that such a correction is sti ll unce rtain 
as discussed by the authors. IBotticella et al.l (|2008r ) did 
a modeling of host gala xy extinction following the recipe 
of lRiello fc PatatI (l2005l ):' their rate measurements take it 

into account both for SNIa and CCSN. 

-3 



We use the SNu or ''SuperNova Unit" , 1 SNu = 1 SN/lO^^LBo/cgntMbdelization of the supernova rate evolution 



2. Supernova rate measurements 

Together with the first high-redshift supernova searches, 
the 1990's have seen the fi rst type la rate measurement 
beyond the local universe ( Pain et al. 19961 ) at z 0.4, 
using for the first time a homogeneous sample of CCD 
discovered SNe. Other supernova searches or surveys have 
provided an estimate of th e rate at va rious redshifts: 
Hardin et all (I2OO0I): iBlanc et al. 



EROS 



Neill et al 
(I2OO7I) 



20041 SCP 



Painet all (l2(302D. S D SS - iMadewick et all (l2003l). HZT 
Tonrv et al.1 (l2003h . iBarris fc Tonrv' f200(h. GOODS - 
Dahlen et all j2004[ l: IKuznc t sova et al. ( .20081). SNLS - 



(I200e 



Subaru Deep Fi e ld - iPoznanski et al 



STRESS - IBotticella et all (l2008l ). Tab. [H sum- 



marizes all the present SNIa rate measurements includ- 
ing the rate in the l ocal u niverse derived by combining 
the ICappellaro et all (|l999f ) determinations for all galaxy 
kinds, for homogeneity with the high redshift measure- 
ments. 

Several attempts have been made to measure the SNIa 
rate in galaxy clusters. This approach is interesting since 



The SN rate (TZsn) as a function of time can be modeled 
as the convolution of the SFH (5*) with the distribution 
function (/sn) of the delay times (r): 



krit - r) • *(i - T)/gN(T)AsN(< - T)dT 



where t is the time elapsed since the beginning of star for- 
mation in the system under analysis; kr is the number of 
stars per unit mass of the stellar generation with age r; 
T'min and Tmax bracket the range of possible delay times; 
AsN is the number fraction of stars from the stellar gener- 
ation born at epoch {t — r) that end up as SN. For core- 
collapse events, the progenitor life-time is very short, (r < 
0.05 Gyr), so that Eq. ([J) can be approximated as: 



7ecc(0 = fcr^cc • *(0 



(2) 



with the SFR and number of SN per unit mass of the parent 
stellar generation (fcr^cc) evaluated at the current epoch. 
For type la events a very wide range of delay times is pre- 
dicted from stellar evolution, with Tmin of a few tens of Myr, 



and Tniax of the order or one Hubble time, or more. We can 
still introduce a little simplification to Eq. ((T|) by assuming 
that the number of SNIa from 1 M0 stellar generation does 
not depend on time and write: 



7^Ia(^) = krAi, 



-)/la(T)dr 



(3) 



which shows that the SNIa rate evolution depends both on 
the SFH and on the distribution of delay times. 

Both ^cc and Ai^ can be evaluated from stellar evolu- 
tion, including a choice for the IMF. However, these theo- 
retical values depend on assumption on the mass range of 
the progenitors. In addition, Aia depends on the distribu- 
tion of binary separations and mass ratios, the outcome of 
the mass exchange phases, and of the final accretion on top 
of the WD. Therefore, we prefer to determine these quan- 
tities empirically from the fit of the SN rates, and compare 
the derived values to the theoretical expectations. 

3.1. The Star Formation History 

The SFH, 'l'(t), is currentl y measured up t o z ~ 6 (with 
some estimates up to 2; 10. iBouwens et ah ( 20051 )) from 
galaxy surveys, using luminosity tracers and exploiting 
cor relations of SFH with UV, H^, FIR luminosities (see 
e.g. lKennicut3(|l998l) ). It has been pointed out that short- 
wavelength measurements lead to an underestimate of the 
SFR because of dust absorption at high-redshift. Then FIR 
luminosity which is less affe cted by the dust obsc uration 



might be a better indicator (jCharv fc Elbadl200lD. Fig. HI 
shows the c ompilation of SFR measur ement by iHopkinsI 
( 2004 ) and Hopkins fc Beacom ( 20061 ). who rcscale the 
literature data to derive a homogeneous sample, with the 
same cosmology, IMF and correction for absorption. Al- 
though the data points generally describe a consistent pic- 
ture, there are large discrepancies at given redshift, which 
likely follow from the different syst ematics in the het- 
erogen eous sample (see discussion in iHopkins fc Beacom 
( 2006[ )). As a result, several laws for '^(t) are consistent 
with the dat a. For our exercise, w e consider the following 
two options: ICharv fc Elba3 (|200l[ ) SFH ~ hereafter CE, a 
model derived from far-IR SFR measurements under the 
form of acceptance interval, which we implement with its 
average values at each redshift ( middle blue curve in Fig.[T]) ; 
and the function introduced bv lCole et all (|200lh and used 
by IHopkins fc BeacomI ( 2006f l to fit their SFR measure- 
ments compilation: p*(2) = 0.7/i7o(a -I- bz)/{l + [z/cY). 
Our fit to the data based on a Salpeter IMF gives: a = 
0.0134, b = 0.175, c = 2.93, d = 3.01. The two SFH mainly 
differ in their peak in redshift, which is ~ 0.8 and 2 re- 
spectively for CE and Cole. They also correspond to two 
different values for the total gas mass cycled into stars, 
with CE being a factor 1.3 larger than Cole's. 

Since ^{t) is derived from sampling the massive stars, a 
large extrapolation factor is implied in the derivation of the 



total SFR, which depends on the assumed IMF. In addi- 
tion, the conversion factor between UV luminosity and star 
formation rate in massive stars depends on the slope of the 
IMF in the high mass range. As long as this slope is simi- 
lar to the Salpeter's the conversion factor does not change, 
and the "^{t) derived for Salpeter IMF can be rescalcd to 
another IMF with a different proportion of low mass stars. 
In our computations we adopt a Salpeter IMF. 

3.2. The Distribution of the delay times 

The distribution function of the delay times describes 
the proportion of early and late explosions past an instan- 
taneous burst of star formation. Indeed, for such a case 
Eq. ([3|) becomes: 



7^Ia(^) = kvAi, -MB-fuir^t) 



(4) 



where A4b is the total gas mass turned into stars in 
the burst. The theoretical /la function is uncertain be- 
cause of various reasons: on the one hand there exist 
different classes of potential SNIa progenitors, on the 
other, the delay time distributions of each class depend 
on the distribution of binary parameters, and on the 
effect on them of the mass exchange phases which oc- 
cur dur i ng th e evolution of close binary systems (see 
Greggiol 20051 for a discussion). In the literature dif- 
ferent formulations can be found for /la, some stem 
ming from stellar evolution prescriptions (e.g. | Sadat et al 



1998 : lYungelsonfc Livioll2000HHan fc Podsiadlowskill2003 



2004 lGreggioll2005h . others from convenient analytical 



pararnetrizations (e.g. 



Madau et al. 1998: Mannucci et al 



2005t Scannapieco fc BildstenT2005l ). others from the fit 



of the observed ti me behaviour of the SNIa rate (e.g. 
Strolger et al.ll2004l ). For this paper we have selected a few 
of these proposed distribution function of the delay times 
to evaluate their impact on the interpretati on of the cos- 



mic ev olution of the SNIa rate, namely the | Madau et ah. 
( 19981 ) formulation, a selection of Greggio (l2005l) models , 



and a selection of the Han fc Podsiadlowskil ( 2003 . 2004) 
simulations. The first are parametric analytic formulae 
which allow the exploration of a wide range of shapes; the 
third is meant to illustrate the result of specific predictions 
from binary population synthesis 



Left panel in Fig. [2] shows examples of the lMadau et al 
( 1998[ ) distributions, i.e. 



/ia(T) OC J exp 

max(mia^^.^,(r/10)-o<') 



T - t„ 
TMada 



-2.35 



dm (5) 



where = lO/rn^'^ is the main sequence lifetime (in 
Gyr) of a star of mass m (in Mq), and TMadau is a charac- 
teristic delay time, considered as a free parameter, mia^j^ 
and mia„,^^ bracket the mass range of SNIa progenitors, for 
which we have adopted 3 and 9 Mq respectively. As TMadau 
increases, the distribution of the delay times becomes wider, 



Q 



still with most of the explosions occurring at early delays. 
For a long TMadau the distribution becomes flat. 

T he centra l pane l in Fig. [2] shows the selected models 
from iGreggio ( 2005 ) . These are analytic delay time distri- 
bution functions based on general arguments on the evo- 
lutionary behaviour of stars in binary systems. The differ- 
ent SD models, shown as dotted lines, correspond to dif- 
ferent schemes for the derivation of the distribution of pri- 
mary and secondary mass in the SNIa progenitors, and to 
different options for the explosion (Chandrasekhar or Sub- 
Chandrasckhar). All cases adopt a flat distribution of the 
mass ratios (i.e. /i = 0). The DD models, shown as solid 
lines, are characterized by different assumptions about the 
result of the first common envelope evolution (WIDE or 
CLOSE systems), and different slopes for the distribution 
of the separations of the DD system at birth. This is de- 
scribed by the /3 parameter, whose definition is different 
in the CLOSE (/3g) and WIDE DD (/3a) schemes. Larger 
values of these parameters imply flatter slopes of the dis- 
tribution of the delay times at late epochs. All DD models 
in Fig. [2] assume that the minimum mass of the secondary 
component in SNIa progenitor system is ~ 2.5 Mq, which 
corresponds to a maximu m nuclear delay of Tn,x = 0.6 Gyr. 
We refer the reader to the lGreggid (j2005l ) paper for a more 
detailed description of the parameters. All of the distribu- 
tion functions in Fig. [2] exhibit an early rise, followed by 
a populated peak of variable width, followed by a tail of 
late explosions. The proportion of early-to-late epoch ex- 
plosions is highly variable among the different models. 

In Fig. [2l we also plot our adopted models from 
Han fc Podsiadlowskil (|2003l l2004h in the right panel. The 
Han DD model (red), refers to double degenerate progen- 
itors, with a standard parametrization for the common 
envelope evolution (acE = ctth = l,aRLOF = 0.5 in 
Han fc Podsiadlowskil |2003f ) notation). The Han SD model 
(green) is a single degenerate, Chandrasekhar model with 
initial secondary masses between 2 and 3.5 M0. These 
limits follow from the requirement of stability during the 
mass transfer from the donor to the WD, so that the latter 
component can reach the Chandrasekhar mass. The direct 
consequence of the lower limit on the secondary mass is 
that this SD model lacks the late delay times tail in the 
/la function. 



4. Results 

In this section we explore which constraints can be de- 
rived by the analysis of the cosmic SN rate, by applying 
Eqs. ^ and Fits of SN rate are done by a method, 
providing an estimate of goodness-of-fit as a "x^ proba- 
bility" . We note that the supernova rate statistics is pois- 
sonian, which converges toward gaussian statistics for big 
numbers. Although most of SN rate used in this paper are 
based on small numbers of events, we assume that we can 
use the fitting method. Systematic errors (when avail- 
able), are taken into account by adding them quadratically 



to the statistical error. For each fit we derive a value for 
the /crAsN factor, i.e. the number of SN per unit masJ^ 
from the parent stellar population needed to account for 
the observed level of t he rates. 

The SNIa rates bv iDahlen et al. ( 2004 ) have been re- 
cently re-computed by iKuznetsova et al.l (j2008[ ) with a 
more sophisticated method. For our fit we adopt the most 
up-to-date values. 

Similarly, f or the STRES S surv ey, we adopt the 
CCSN rate in iBotticella et all (|2008h which updates the 
ICappellaro et all (|2005h " value on the basis of a larger SN 
sample. 

4.1. SNIa rate and SNIa progenitors 

Fig. [3] illustrates our fit with the Madau's functions. The 
best case is obtained with TMadau = 1 Gyr, corresponding 
to a very peaked distribution of the delay t i mes (c f. Fig. ^ , 
in combination with CE SFH. Cole et all (2001)'s cosmic 
SFR is excluded by about 3-sigma (V{x^) = 5 • 10 "^). Sim- 
i larly, Fig. m shows the fit to th e data with the Greggiol 
120051 ) distribution functions. The lCole et aD (l200l|) 's SFH 
do not provide a very good fit, while they all get a signifi- 
cant fit probability (over a few percents) when CE's SFH is 
adopted. Fig. [5| shows the fit based on our Han and Podsi- 
adlowski selected models. Again the best fit is obtained in 
combination with CE SFH. The main difference between 
the two SFH is the p osition of maximu m, Zma.y ^ 0.8 for 



Charv fc ElbazI (|200l[ ). and z^ax ~ 2 for lCole et al.l (|200lh . 

We conclude that the SNIa rate data favor a maximum of 
S FH at late epoch in t he history of the universe. 

F5rster et al. ( 20061 ) obtained a similar conclusion with 



different ingredients: the typical delay time of SNIa pro- 
genitors inferred from the cosmic evolution of their rate 
depends on the SFH adopted. Inspection of Fig. 21 how- 
ever, shows that even for a given cosmic SFH, the model for 
the SNIa progenitors is not well constrained, with all mod- 
els running very close to each other. Although an adopted 
cosmic SFH allows to solve for one most probable delay 
time, this does not tell much about the progenitors of SNIa, 
whose delay times are spread over a wide distribution. 

In general, the diagnostic power of the cosmic SNIa rate 
is relatively weak, due to the fact that these data are aver- 
ages of the SNIa rate in all galaxy types, over the SFH of the 
whole universe. In other words, it results from the convolu- 
tion of the fia function with a wide ^'(t) function. A better 
way to constrain the SNIa model would consist in consid- 
ering the SNIa rate evolution as function of galaxy type or 
in galaxy clus ters, which are characterized by a relatively 
narrow "^{t) (jCal-Yam et al. 2002). This is attempted in 
the next section. 

The values obtained for the fcr^ia factors (where fcr = 
2.83 Mq^ for a Salpcter IMF) from fitting the cosmic SNIa 



^ This is the mass of the stellar population at zero age, which is 
larger than its mass at any subsequent age because of the mass 
return from stars. 



A 



rate with Greggio's distributions range in [0.45 — 0.55] 



10 



-3 and [0.58 - 



0.71] • 10-3 Mq^ when adopting CE 
and Cole SFH respectively. The difference reflects the dif- 
ferent total star formation for the two options, with CE 
providing more stellar mass than Cole's, so that the level of 
the SNIa rate is recovered in combination with a lower real- 
ization probability of the event. The fit obtained with Han 
distributions and CE SFH yields similar values: fcr^ia = 
[0.33; 0.53] • lO'^ Mq^ for DD and SD models respectively. 
Notice that the fits obtained with Han's distributions in 
Fig. [5] bracket the curves obtained with Greggio's distribu- 
tions. We compare these derived factors to the number of 
stars per unit mass of the parent stellar population with 
mass between 3 and 9 Mq , which is an indicative range of 
the mass of SNIa progenitors (Umeda ct al. 1999). Since 
Salpeter IMF provides 0.022 Mg^ such stars, the fcr^ia fac- 
tors require that about 2-3 % of them give rise to a SNIa 
event. 

We notice the lar g e disc repancy between rates measured 
by iBarris fc To"^ (|2006h at z > 0.5 and by othe rs: at 
z 0.55 the discrepancy with the Pain et ^dl ( 2002 ) mea- 



surement reaches more than 3cr. Since the lBarris fc Tonrv 



( 2006f ) measurements are not fully based on spectroscopi- 
cal confirmation and do not include systematics errors, we 
also perform the fit excluding the 3 discrepant points. The 
resulting fcr^ia factors change for less than 20 % but the 
quality of fit drastically increases, 91 to 97 % for Grcggio 
models combined with CE SFH, and 54 to 70 % for Greggio 
with Cole SFH. This favors Han's DD model (97 % and 74 
% for CE and Cole SFH respectively) , while their SD model 
gets a 50 % fit probability. The quality of the fits has become 
so close that we cannot conclude for one or the other SFH 



model . In this respect, we notice th at rate bvlDahlen et al 
(|2004f) at z ~ 0.8 is similar to t he iBarris fc Tonryl hOoT 



value, but its redetermination bv'Kuzne tsova et al.l 
and the new data bv [Poznanski ct al.i (2007) favor lower 
SNIa rates at these redshifts. 

4.2. SNIa rate in galaxy clusters 

The galaxy population in clusters is dominated by Early 
Type Galaxies (ETGs), which, accordin g to cu rrent under- 
standing, contain mostly old stars (see iRenzin i 2006, for 
a review). This is so not o nly in the local univ erse, but at 



least all the way to z 1 ( Tanaka et al. 20051 ) . Therefore 



it is generally expected that the evolution of the SNIa rate 
in clusters of galaxies basically traces the evolution of the 
SNIa rate in ETGs, which closely scales as the distribution 
function of the delay times, as l ong as their SFH can b e 
described as a short initial burst (|Gal-Yam fc Maozll2004l ) . 

In a cluster, the SNIa rate in SNu, 7?.ia,sNu, is given by 
the sum over the cluster members of their individual rates 
in SNu, i?ia,SNu, weighted by their contribution to the total 
B luminosity, fs- 



"^la.SNu — ^ ^Ia,SNu,i ' fs.i- 



Both the individual rates and fs^i factors vary with red- 
shift, so that in order to compute the expected trend for 
the rate in clusters one should specify the galaxy popu- 
lation in detail. On the other hand, in spite of the ob- 
serve d variation of the cluster galaxy pop ulation with red- 
shift (jDressler et al.lll997trPoggiantii2006h . the evolution of 
the bright portion of the luminosity function in clusters is 
consistent with that of a passively evolving model in which 
galaxies form their stars in a sin gle burst at high (z ~ 2, 3) 
redshift ( de Propris et aLlfiggoh . This justifies approximat- 
ing the right hand side of Eq. ([6]) with the evolution of 
the rate obtained for such stellar population. By applying 
Eq. ([4]) we have: 



7^Ia,SNu(^) ^ krAu ■ (M/Lb) ■ fuir ^ t - tf) 



(7) 



where t is age of the universe at a given redshift, and tj is the 
formation epoch of the stellar population. Thus, the SNIa 
rate in galaxy clusters as a function of redshift (in SNu) 
scales as the product of the mass-to-light ratio and the dis- 
tribution function of the delay times, both of which depend 
on the age of the stellar population at the various redshifts. 
The fit of this relation to the corresponding d ata yields a 
value for the fcr^ia factor. ISharon et al.l (|2007l ) present the 
SNIa rate in clusters as a function of redshift. We plot their 
data in Fig.[6l together with our fit obtained with a for- 
mation redshift of 3. We checked the dependence of the fit 
on formation redshift by adopting Zf = 2,4 and found it 
to be small. We have included in the fit the local value as 
measured bv lCappellaro et al ( 19991 ) for E/SO galaxies in 
the field. Actually, even if ellipticals in th e field may be 



slight ly younger than those in clusters (e.g. iThomas et al 



200i), the low sensitivity of our fit to the formation red- 
shift makes th e use of this local meas urement reliable. If 
we exclude the Cappellaro et al. ( 1999h point from the fit, 
we derive a ~ 40 % larger value for both /cr^ia and its 
error bar. The mass-to- light ratios as function of re dshift 
have been derived by integrating iGirardi et al. ( 2002h solar 
metallicity isochrones with a straight Salpeter IMF. 

Our fit with Madau's functions shows that the data 
favour relatively large values of TMadau^ the fit with 
'Hviadau = 1 implies a too rapid growth of the SNIa 
rate with redshift . This result is different from that in 
Maoz fc Gal-YamI ( 2004[ ). who, based on the same data, 
conclude for a relatively short TMadau (< 2 Gyr for a forma- 
tion redshift of 2). The different conclusion likely follows 
from the different approach; in particular we consider the 
variation of the mass to light ratio with redshift, while they 
do not. The decrease of M/Lb with increasing redshift 
counterbalances the growth of the distribution function of 
the delay times, resulting into a slower expected evolution 
of the rate in SNu with redshift. 

Fig. [6] shows that due to the large error bars, no real 
constraint on the distribution function of the delay times 



(6) 



^ The fit for Han SD model cannot be performed, since the SNIa 
rate from this model is zero when the population is older than 1.5 
Gyr (see Fig.|2}. 



can be derived at prese nt. Most of the error bar reflects the 



low number statistics ( Sharon et al.l 120071 ) so that future 



surveys will greatly improve the situation. The models in 
Fig. [6] illustrate the diagnostic from this kind of data: if the 
rapid growth of the SNIa rate with rcdshift is confirmed, the 
DD models in Greggio's formulation would be unfavoured, 
since they correspond to a flat behaviour of the rate at red- 
shifts smaller than ^ 1. The steep trend exhibited by the 
SD Chandrasekhar models (blue and green dashed lines) 
mirrors the steep late epoch decline of the distribution func- 
tion of the delay times for these kind of progenitors, due to 
the requirement that the WD reaches the Chandrasekhar 
mass by accreting the e nvelope of a pro gressively less mas- 
sive companion (see e.g. lGreggidl2005l for more details). 

The derived values for the number of SNIa from 1 Mq 
stellar population are affected by a larger uncertainty with 
respect to what we get from fitting the rate in the field. 
Averaging over all Greggio's models we get fcr^ia = 1.7 • 
10-3 MQ\notfarfr om those obtained in the previous sec- 
tion. We acknowledge a difference of a factor of 3, favour- 
ing events in clusters, at a ~ 2 cr level. Given the different 
ingredients used in the two computations, assumptions on 
the global SFH and age distribution, and the uncertainty 
on the observed rates, we regard this result as encouraging, 
although the discrepancy needs further investigation. 

An independent constraint on the /cr^ia factor can be de- 
rived by considering the total Fe budget of galaxy clusters. 
Assuming that every SNIa event ejects an average mass mpc 
of Fe, their contribution to the Fe Mass-to-Light ratio is: 



Fcja 



= TOR 



(8) 



where M/Lb is the ratio between the total gas mass turned 
into stars and the current B band luminosity of the typi- 
cal elliptical in the cluster. For a Salpeter IMF, M/Lb — 
13.5 Mq • Lq^ at an age of 11.4 Gyr, corresponding to a for- 
mation redshift of 3. The total Fe mass to light ratio in clus- 
ters is ^ 0.02 fe.g. lDe Grandi et al.|[200l : ~ half of which 
is contributed from SNIa (see Sect. 5, also Greggio 2008, 
in pr eparation). Adopting mpc = 0.75 Mq (jlwamoto et al. 
I999I summing the ejected masses of ^^Fe and ^^Fe), we 



get ferula ^ 10" Mq . We conclude that the three inde- 
pendent estimates of the A;rj4ia agree within a factor of 2. 
The uncertainties currently affecting each of them prevents 
a more accurate determination. 



4.3. CCSN progenitor masses from CCSN rate 

Due to the short lifetime of their progenitors, in a stel- 
lar system the CCSN rate is directly proportional to the 
current SFR (Eq.[2]), with: 



kvA 



mcc„ 



r^cc 



r 



m^{m)dm 



(9) 



where $ is the IMF by number, [m*^;^, m*^^^] is the whole 
stellar mass range and mcc„i„ and TOcc„ax are respec- 



tively the minimum and maximum stellar mass leading to 
a CCSN. 

For a Salpeter IMF between 0.1 and 120 M©: 



fcr^cc = 0.126 • (m 



-1.35 ™-1.35 



CC„ 



CC„ 



(10) 



By fitting Eq. ^ to the redshift trend of the CCSN rate 
with the two SFH laws described in Sec. 13. II we obtain the 
related values of the constant fcr^cc- Although the scaling 
of the SFH to the CCSN rate depends on the IMF, we only 
consider Salpeter IMF because this has been adopted for 
the derivation of the SFH laws. The fitting values of kpAcc 
are reported in Fig. [7l giving similar values for CE and 
Cole SFH. For a given SFH we get a bijective relationship 
between rncc^i„ and mcc„^^^-, which is very sensitive to 
fcr^cCi as plotted in Fig.[8l 

The minimum mass TicCmi,, is actually set by the heavi- 
est star that produces a wh ite dwarf, which ca n vary from 
6 Mq to 11 Mq according to lHeger et"aI1 (|2003h . They used 
9 Mq as an average value. It is supported by the obser- 
vations of SNII events: the lowest SNII progenitor mass 
known so far has been derived from the type II-P event 
SN 2 003gd as Stj Mq (jSmartt et al.ll2004l : Hendrv et al 



200i), though careful modeling of the observations favors 
a progenitor mass c loser to the upper limit of this range 
(jHendrv et al.ll2005[ ) . The maximum mass TOcc,„ax much 
less constrained, but should be around ~ 40 Mq, or less. 
The highest main sequence masses that have been derived 
for CCSN event arc for some type Ic connected with a hy- 
pernova event, with TOcc,„ax ^ 40 Mq. This kind of SN 
is of course included in the generic class of "CCSN event" 
for which the rate history has been measured. For instance 
the type Ic SN 20031w has a progenitor with a mass within 
[40 - 50 Mq] (jMazzafi et al. 2006). Within the type II-P 
SNe class, SN 2003Z has t he highest "observed" mass with 
mcc„,x = 36.5 ± 4.9 Mq (|Pastorelloll2"oo5lZampieri et al. 

Fig. [8] shows that the observed and theoretical rncCmin 
and JTicCmax are only compatible with the high side of the 
1-cr error interval of krAcc for both CE and Cole SFH. 
The favored ?7T,cc„i„ is TO.cc,„i„ ^ 10 Mq to account for the 
observed wcc^ax ■ since CCSNe explode in star forming 
regions, a fraction of the events is likely to escape detection, 
so that the rate in Fig. [7| would be a lower limit to the 
real CCSN rate (see sec. [2]). In this case, the fit with the 
SFH would lead to a larger value for krAcc, which could 
accommodate a wider range of CCSN progenitor's masses. 



5. The released iron density evolution 

The iron in the Universe is produced or released in the 
ISM by supernovae. Knowing the supernova rate history 
and the distribution of produced iron, we can estimate the 
evolution of the iron density with redshift. As a function of 
cosmic time, it is generally given by: 



PFc(i) = / MFcit')-nsN{t')dt' 



(11) 



where A4fo is the distribution function of iron mass released 
per event and tic is the epoch of rcionization corresponding 
to the first stars, and then the first SNe, and first iron 
released in the ISM. 



5.1. From type la SNe 

A SNIa releases in the interstellar medium (ISM), on av- 
erage, 0.5 Mq of ^^Ni (that becomes of ^^Fe after radioac- 
tive decay). During the explosion others isotopes of iron 
are produced, ^^Fe, ^^Fe and ^^Fe. The present approach 
is sensitive only to the "observable" iron which powers the 
light curve of SNIa, that is ^^Fe. Then we expect to miss 
10 % to 15 % of iro n according to delaye d detonation nucle- 
osynthesis models ( Iwamoto et al.|[l999() : for 3-dimensional 
pure deflagration models, the proportion of iron isotopes 



other than ^^Fe varies from 15 % to 25 % ([Travaglio et al 



To take into account all iron, we assume that ^^Fe 
represents 85 % of all iron released by SNIa in the ISM. 
Then, in Eq. (jlip the iron distribution is given by Adpe = 
1.18 A^56p(, = 1.18 TWseNi, assuming that this ratio does 
not evolves with cosmic timj^. 

The distribution of released ^^Ni m ass is taken from 
Leibundgut (|200Gl ). lMazzah et al.l(|200Gl ) and lStritzinger et al 
(l200d) : it is plotted on Fig.[51 The averaged ^^Ni mass as de- 
rived from observations, 0.5 Mq per event, l ies between the 



predic tion of DD models (0.6 M© to 0.8 M^ jlwamoto et al 



( 19991) ) and predictions of pu re 3-dimensions deflag ration 
models (0.24 M© to 0.44 Mrr. . ITravaglio et all (|2004[ )). 

To render the SNIa rate measurement and integrate it 
easily, we fit it by an exponential between z = and 0.59, 
and by a straight line for z > 0.59 as shown on Fig. (TUl 
Eq. (fTTj) is integrated by a Monte-Carlo method in order to 
keep the observed shape of the nickel distribution, under 
the assumption that the distribution of ^^Ni as shown on 
Fig. [9] is representative and independent of cosmic time. 
Beside the ad hoc shape for the SNIa rate, we use the best 
model fitting the rate evolution, that is the Greggio SD- 
SCh distribution convolved with the CE SFH (Fig.|4|). The 
result is plotted on Fig. [TOl 

Since the model produces SNe la at redshift larger than 
2, for this case we get a larger amount of integrated 
iron, with respect to the ad hoc fit. Still, the difference only 
amounts to a factor of 1.2 at the present epoch. So SNe la 
in the distant universe have little influence on the total iron 
produced. 



^ Actually this rat io probably evolves with redshift since 
ITravaglio et al ] l|2005h found an anticorrelation of A^sepg/Afpo with 
metallicity, going from 0.93 at 0.1 Zq, to 0.83 at Zq and 0.70 at 
3 Zq . But such a study is beyond the scope of this paper. 



5.2. From CCSNe 

Like for SNIa, only powers the light curve of CCSN, 
by rad ioactive decay to ^^Fe. According to llwamoto et al 



1999f) a typical CCSN ejects 93 % of the total iron as 
^^Fe; thus in our total computation, we take Aipe = 
1.08 M56Fe = 1. 08 M56N\ in Eg . (fTTI) 



According to IZampieri et al.l ([20031) and iHamuvl (|2003l ) 
the ^^Ni mass released by SNII is distributed as shown in 
Fig. H based on 29 SNe II. The average Ni mass is 0.066 
Mq . This distribution covers more than two orders of mag- 
nitudes, while the distribution of Ni produced by SNIa 
spreads less than one order of magnitude. Eq. pT|) is inte- 
grated over this distribution by a Mont e-Carlo method, set- 
ting t he redshift of the first stars to 20 (jFuUer fc Couchman 
l2000f) . The CCSN rate is modeled by two SFH models as 
previously discussed. They are plotted on Fig. [Til together 
with the results of the integration. The two SFH models, 
once calibrated on the observed rates, produce very close 
Fe production curves. 



5.3. Total 

Fig. [12] shows the total density of iron produced, sum 
of the contribution of SNIa and of CCSN. It is interesting 
to note that today's iron density comes from half-part of 
SNIa (Fig. [ini) and half-part of CCSN (Fig. [II]). Indeed, 
on average one CCSN releases 10 times less Fe than one 
SNIa. However, the CCSN production per unit gas mass 
transformed into stars is about 10 times that of SNIa pro- 
duction, as result from the ratio fcrAcc/^r^ia ~ 7.5. We 
remark that the evolution of the ir on density does not 
depen d strongly on the SFH model. ICalura fc Matteucci 
( 2006h present models of cosmic metal production based on 
chemo-photometric models of the evolution of dwarf galax- 
ies and spheroids in the universe. The models detail the evo- 
lution of various chemical species in different phases (ISM, 
IGM, Stars), and, in particular, that of the total Fe den- 
sity. Fig. [T2l shows as blue lines the predictions of Calura & 
Matteucci models for two kinds of galaxies, irregulars and 
spheroids, and their sum. Since the contribution of disks is 
missing, we cannot compare directly the Calura and Mat- 
teucci predic tions (solid blue line) with our empirical curve. 
However, in ICalura fc Matteucci ((20041 the iron produc- 
tion rate by disks becomes important only at relatively low 
redshift {z < 1.2), and at zero redshift the contribution of 
disks is ~ 0.6 that of ellipticals. W e may then extrapolate 
the results of ICalura fc Matteucci (l2006l) into a total ex- 
pected Fe production in the Universe of 8 • 10^ M© Mpc~'^. 
This is in excellent agreement with our estimate, which is 
Pfc = (8.5 ± 0.9) • 10^ M© Mpc~^, as averaged from our 
various estimates. We also expect that the corrected curve 
for the total Fe in Calura and Matteucci's models is similar 
to what plotted in Fig. [12] for redshift in excess of ~ 1.2, 
so that also in the models a large fraction of Fc in the uni- 
verse is produced at z < 1. Given the completely differ- 



ent approaches for determining the cosmic increase of total 
Fe density, we consider the models and our empirical es- 
timate in very good agreement. Our total current Fe den- 
sity is also in excellent agreement with the value of ppo = 
8.6 • 10^ Mq Mpc"^ bv lFukugita fc Peebles! (|2004l ). 

We label the right axis in Fig. [12] with the Fe abun- 
dance relative to solar, having considered Zpc = Pfc/pb, 
where pe = Pc, ■ i^B, Po = 1.36 • 10" /i^gMoMpc"^ is the 
critical density of the universe and 1 7^ is the density o f 



baryon (f^s = 0.049 ± 0.002 /ifff, see ISpergel et al.ll2003l ). 



which gives: Zpo = (1 .28 ± 0.14) ■ 10"^. Adopti ng the so- 
lar Fe abundance by iGrevesse fc Sauval (1998) ^fg,© = 
(1.3 ± 0.1) ■ 10"'^, we obtain an average Fe abundance of 
the universe at zero redshift of Zpc = (0.10 ± 0.01) Zpc,Q, 
which is a factor of ~ 3 -4 lower than the Fe abundance in 
clusters of galaxies (e.g. lRenzinll2004l ). 

Does this mean that the stellar production of Fe is dif- 
ferent in clusters from that in the field? To answer this 
question we consider the ratio between the Fe and the stel- 
lar density 77F0 = PFo/p*- In clusters: rjFc.ci = ^FcJCM • 
MfC M I -I- Zpn. * — 3.3 ^Fo,o having adopted M,,/Mtcm = 
0.13 |Ettorill2007h . Zfc,icm = O.3Z0 (jXamura et al.ll2004h . 
and a solar iron abundance for the stellar component. By in- 
tegrating our SFH and correcting for the mass return from 
SNe and stellar wind£l], we find a current ~ (8.2 ± 4.1) • 
10^ /17OM0 ■ Mpc"^ which gives 77F0 = (0.8 ± 0.4) ■ Zfo,© 
a factor ~ 4 lower than in clusters. However, while these 
estimates of 77fc,c1 and of the Fe density in the field do not 
depend on the IMF, the stellar density in the field does. 
A more realistic IMF, with a turn-over at the low mass 
end, implies a lower stellar density for the same amount 
of UV photons. For example, the mass fraction in stars 
more massive than 1 M© is 0.61 and 0.39 for Chabrier 



( Chabrieiil2005f ) and Salpeter IMFs respectively; the mass 
fraction returned to the medium from stellar winds and 
SNe is - 0.45 and 0.28 for the two IMFs. Therefore, for 
a Chabrier IMF, the current stellar mass density becomes 
^Cha = (4 ± 2) ■ 10^ MoMpc"^ to yield r^p^cha = (1.6 ± 
0.8) ^Fe,0, close to the value in clusters. In addition, the 
rate of CCSN may be affected by a more severe incom- 
pleteness than adopted in the values in Tab. [21 particularly 
at high redshift, whe re some events could occur in heav- 
ily obscured regions ( Mannucci et al. 2003[ ). We conclude 
that the Fe production from stars in the field and in clus- 
ters can be brought into agreement, if the IMF has a turn 
over at the low mass end. The different abundance, then, 
stems from a different efficiency of baryon conversion into 
stars p'x^^ I Pb^ which adopting the figures above, is ~ 0.12 
in clusters and ^ 0.06 in the field. 



6. Conclusions 

In this paper we have investigated on the constraints 
which can be derived from the cosmic trend of the SN rates 



on the SFH, on the distribution of the delay times of SNIa 
progenitors and on the mass range on the CCSN progen- 
itors. In addition, we have determined the local Fe abun- 
dance as resulting from the integration of the observed SN 
rates. Our results can be summarized as follows. 

The observed cosmic evolution of the SNIa rate has been 
fit by adopting two laws for the SFH which bracket mo st of 
the observational data, na melv lCharv fc Elbaj (|200l[ ) and 



Cole et al 



^ for a Salpete r IMF the correction factor amounts to 0.72 (e.g. 
ICole et al.ll200ll) 



( 200l[ ). in combin ation with various d istribution 

functions o f the delay time s: iMadau et al. (1998) models, a 
selecti on of Greggiol (12005) ni o dels, a nd a s election of mod- 
els by iHan fc Podsiadlowskil (|2003l |2004| ). We find that, 
while, within the error, it is hard to distinguish among the 
possible progenitor models , the fits are fair ly sensitive to 
the adopted SFH (see also Oda et al. l2008h . The increase 
of the SNIa rate with redshift is better reproduced with 
C E rather than Cole SF H. However if rate measurements 
bv lBarris fc Tonrvl (|2006l) were systematically too high, the 
two SFH would be rather equivalent in fitting the SNIa rate 
data. The inability of the cosmic SNIa rate to constrain 
the progenitor model is related to its being an average over 
a galaxy population which spans a wide range of proper- 
ties. Measurements of the SNIa rate in different kinds of 
galaxies would help us to relate the observed rate, and its 
redshift evolution, to the parent stellar population, in turn 
providing better constraints on the shape of the distribu- 
tion function of the delay times. 

As an attempt in this direction, we have considered the 
redshift evolution of the SNIa rate in galaxy clusters, as- 
sumed to harbor mainly old stellar populations. Due to 
the large error bars affecting the current measurements, all 
SNIa progenitor models provide acceptable fits. Madau's 
models with TMadau shorter than 1 Gyr are however 
excluded. Among Greggio's models, the SD family seems 
favoured, but we emphasize that the error bars on the data 
are currently too large to draw any conclusion. 

The fit to the cosmic SNIa rate and to the rate in clusters 
provides the number of la events per unit mass of the parent 
stellar population (fcr^ia) which is needed to explain the 
observed level of the rate. The field cosmic rate and the 
rate in clusters indicate similar (albeit non equal) values 
for this constant. On the average we conclude that in both 
field and cluster galaxies we have ^ 1 SNIa explosion every 
1000 Mq of stars in the parent population, corresponding 
to ~ a 5 % probability of the SNIa channel, if the primaries 
in the binary progenitors have masses in the range 3 to 9 
Mq (for a Salpeter IMF). 

The redshift evolution of the rate of CCSN is compati- 
ble with both CE and Cole SFH. The favored progenitor 
minimum mass is > 10 Mq, even if some incompleteness in 
the observed CCSN rate could accommodate a (somewhat) 
lower CCSN's progenitor masses. 

The convolution of the observed cosmic SN rates with 
their ejected Fe mass allows us to compute the Fe den- 
sity of the universe at zero redshift. We find pFc = 
± 0.9) ■ lO'^ Mq Mpc~'^, in v ery good ag reement 



with iFukugita fc PeeblesI (|2004l ) and ICalura fc MatteuccH 



The corresponding Fe abundance in the local uni- 
verse is 0.1 ZFe,Q, which is a factor of 3 lower than the Fe 
abundance in massive clusters. We show that, if the IMF 
has a turn over at the low mass end, this difference does 
not imply a different efficiency of Fe production from stars 
in the two environments, but rather it reflects a different 
efficiency of baryon conversion to stars. 
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Fig. 2. Delay time distributions used: Madau ct al. (left), Greggio (middle) and Han and Podsiadlowski (right). They arc all normalized over 
the range 0-12 Gyr. 



1 o 



^IQ-s Chary & Elbaz's SFH and Madau's fla 



o 

Q. 



o 

CO h- 



SFH Chary & Elbaz: 

^Madau = 1 (0-04) 

SFH Chary & Elbaz: 
W = 2Gyr(0.01) 

SFH Cole et al: 

= 1.5Gyr (0.005) 




Redshift 
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Fig. 6. SNIa rate measurements in clusters fitted with Han's and Madau's models (left), and Greggio's delay time dis tributions for an epoch 
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Fig. 10. Density of iron released by SNe la (in red, right scale) for two models of SN rate evolution. The SNIa rate measurements are 
superimposed (left scale) and fitted by an exponential between z = and 0.59, and by a straight line for z > 0.59, or modeled by the Greggio 
SDsubCh delay time distribution convolved with the CE SFH (in black, left scale). 



on 



x10^ 




Redshift 

Fig. 11. Density of iron released by CCSNe (in red, right scale) as computed using different SFR models to fit the CCSN rate (in black, left 
scale). The first stars formation redshift is set to 20. 
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Published rcstframc type la supernova explosion rate measurements. These rates are given together with the mean redshift of the observed SNe 
and the cosmological model assumed (especially for distant values). The fifth to eighth columns give respectively: the number of supernovae 
from which the rate is computed; the effective surface of the survey, or the number of galaxies surveyed; the limiting magnitude of the survey 
and the fraction of spectroscopically confirmed supernovae. NOTES: a) Quoted error bar are first statistics, and second systematics; where only 
one is present, the systematic is not available; b) the value per comoving volume is derived from the value in SNu times the 2dF luminosity 
density (the error on the luminosity density is added q uad ratically to the systern atic uncertainty if distinct); c) For a (0.3, 0.7) cosmology 
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(2005) luminosity function to trans late their rate in SN u; g) statistical and systematic errors are 
i 2008h include the SN sample from lDahlen et al.l ||2004|) ;o compute their rate. 
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Tabic 2 

Published restframe CCSN explosion rate measurements. These rates are given together with the mean redshift or redshift range of the 
observed SNe and the cosmological model assumed. The fifth to eighth columns give respectively: the number of supernovae from which 
the rate is computed; the effective surface of the survey, or the number of galaxies surveyed; the limiting magnitude of the survey and the 
fraction of sp ectroscopically confirmed supernovae. Notes: a) The rate is translated into volumetric units b y using a luminosity d ensity of 



. , o y . 

galaxies as in ICappellaro et aT| l|2005 ) b) W e give here the value corrected for extinction by the authors c) iBotticella et al] l|2008l ) include 
the SN sample from lCa'ppeUaroet'^ 1 1|2005|') to compute their rate. 



